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ABSTRACT 


The hydrodynamical-numerical model of Walter Hansen is 
used to compute tidal heights and tidal currents in the 
Gulf of Maine. The model uses two adjacent open boundaries 
Me which the tides are prescribed at’ cach tame step, using 
jmour tidal constituents. The erid, size ds six nauueecal 
miles and the time step is thirty-one seconds. Seven data 
runs are reported; one uses the tides and no wind and the 
remaining six use uniform wind fields in consort with the 
tides. <A modified method of handling the topographic data 
is used. The pertinent results of the study are: (1) the 
use of Hansen's Model with adjacent open boundaries produces 
broad subjective agreement with observed data, (2) the 
modified method of handling topographic input data is 
workable, (3) wind direction and velocity produce slight 
mations in tidal height, (4) wind direction and vellocivy 
modify the direction of the tidal currents considerably a 
and produce some elo an ite erie ery eae alec oi eal) ee cetach cee current 
speed, and (5) the modifying influences of wind Pledae on 
mrdal heights, tidal current velocities, and tidal ycunwene 
G@irections are more noticeable in shallow water areas than 


in regions of deep water. 
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IT. INTRODUCTION 


The knowledge of the currents of the ocean, in 
adjacent and marginal seas, and in estuaries has long 
been of the utmost importance to navigators and inhabitants 
of coastal regions. Many attempts have been made to under- 
stand the currents produced by physical processes such as 
storm surges, tides, and ocean currents. Years of 
observations and measurements have resulted in a vastly 
increased knowledge of the physical processes of the sea. 
This knowledge has been applied in navigation, coastal 
engineering, commercial fishing, and in prediction of 
tides and tidal currents. There is a need for more 
thorough understanding of the physical processes of the 
ocean in order to cope with such modern problems as 
pollution, deep water salvage operations, mining the 
continental shelf, and farming the sea. 

The need for this knowledge leads to questions 
concerning the possibilities of reproducing oceanic motions 
and how best to accomplish this goal. <A successful model 
for reproducing oceanic =e Ae leads to a capability to 
predict these motions and hence, hastens the possibility 
for various groups of people to benefit from the sea. 

Oceanic motions produced in a numerical model must 
be compared with actual measurements. If the correlation 
is reasonably close, then soni Gonltidence may berplaced 
tim the model. “Atter calibration, the model can be used 


to forecast with all the practical applications that this 


‘> 








implies. Electronic computers make the use of numerical 
models practical and reasonably economical. Numerical 
models are flexible in that the boundary conditions and 
variables involved are easily altered to meet changing 
needs. Hence, more extensive use is being made of 
numerical models in predicting physical motions of the 
sea. 

The ea ody Haiibe altho reer anieanl equations, ina 
general form, provide a background for understanding the 
problems of oceanic motions. These equations form the 
basis for the model developed by Dr. Walter Hansen at 
the Institut fur Meereskunde der Universitat Hamburg 
im 1935. 

Hansen's model has been used to reproduce physical 
processes in a number of bays, tidal estuaries, and 
marginal seas. Examples include the North Sea, Baltic 
sea, The English Channel, The Straits of Gibraltar, The 
Bosphorus, Persian Gulf, Chesapeake Bay, Long Island 
Sound , eames Hawwor., Lonkin Gulis and the Gulliver 
ot. Lawrence. Practically speaking, there are a large 
number of ocean areas where there is a need to determine 
oceanic physical processes which could be solved by 
numerical-hydrodynamic means. 

The primary purpose of this study is to extend the 
coverage of the Hansen model to the region of the Gulf 
of Maine and to verify the model in this region. The 


determination of the general tidal circulation in this 
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region and the modifying effects of various wind fields 
on the tides and tidal currents is an objective of major 
importance. This is true for locations where data for 
comparison are available as well as for points such as 
those on the continental shelf and slope where no such 
data could be obtained. 

The region of study includes the area between 40° 
amd 45° North latitude and extends from Go°30" ta 7 1290n 
West longitude, a total area of 90,000 nautical miles. 
(Figure 1) This primarily oceanic area extends from the 
southern extremity of the Bay of Fundy near Grand Manan 
Island off the coast of Maine south to include Cape Cod, 
Nantucket Shoals, and nearby environs. 

The bottom topography is irregular in the central 
Wilkinson Basin region. (100 fathom contours shown in 
Figure 1) It becomes more regular to the north approaching 
the Bay of Fundy, to the south in the Nantucket Shoals 
area, and as the coastline 1s approached. The underwater 
landscape is principally continental shelf, but in the 
lower southeast corner, the continental slope is crossed 
with an accompanying plunge in water depths from 100 
fathoms or less to 1600 fathoms and deeper. 

The land-sea border is rugged and frequently indented. 
Siais is particularly true along the coast of Maine. The 
various bays and inlets form an interesting and unusual 
coastline which includes several unique features. The 


most dominant feature is that of Cape Cod which extends 
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like a giant fishhook into the sea. The canal across 

Cape Cod is another unusual feature of the area as the 
tides and tidal currents at the eastern and western ends 

of that body make an interesting comparison. several large 
islands such as Martha's Vinyard, Nantucket, and Grand 
Manan are all strategically placed to produce an influence 
upon oceanic motions that may be of interest. Finally, 

the large ee constituents of the Bay of Fundy affect 

the Gulf of Maine region. 

The Gulf of Maine is important for Se eel reasons. 
There are large commercial fishing interests in the area. 
ine once of Boston, Massachusetts, is an important 
Maritime region as is the southern portion of the gulf 
where sea lanes carry trade between New York City and 
various world ports. Numerous hazards to maritime safety 
are attested to by the number of shipwrecks in the area, 
especially on Nantucket Shoals. There are many other 
reasons for the region's importance but these are perhaps 


the most important. 
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ik. THE EQUATION: 


The classical problem of the tides in the world 
ocean can be stated: Determine the tides and tidal 
currents from hydrodynamical equations using the known 
tide generating forces and the geometry of the sea bottom 
and coastline. Tt is then possible to compute the tides 
and tidal currents quentata bagel, using the hydrody- 
Namical differential equations without resorting to 
measurements of these quantities. 

This problem has been studied since the time of 
Newton. Mathematicians and physicists such as Laplace, 
Bermoulli, Poincare, and others, have worked on the tidal 
problem. The common goal of their combined efforts was 
to solve the hydrodynamical equations analytically. 
solutions were obtained only for geometrically simple 
oceans. Recently, Proudman and Doodson developed a 
method of computing ocean tides for regions bounded by 
meridians. In all these efforts, the essential influence 
of the depth distribution on the tides and tidal currents 
was considered in a geometrically simple manner. Con- 
Pocerangle diiferences have been eatnd to exist between 
iepaesS and tidal currents in the actual ocean and those 
computed for geometrically simple oceans. 

The equations used in Hansen's model are based on the 
hydrodynamical differential equations derived from the 


Navier-Stokes equation. They are the equations of Mor LO: 
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and the conservation of mass equation: 


$342 i (Ha) +5 (Hv) = O (3) 


where xX and y are the space co-ordinates; t is the time; 
u and v are the components of the mean velocity; H is 
the total depth; <a water elevation in relation to the 
un-disturbed surface of yy sea; X and Y are the external 
body forces; qe) i Cy) the wind stress components; 
g the acceleration of gravity; f the coriolis parameter; 
xr the coefficient of friction between water and the sea 
imEOOr ; Ay the ,Coetriacient vol “horizontal eddy difttuwsay26.,. 
and fy the Laplace operator. 

Tie sedayvatithusivity vemmis sarc Inver pre tededemats 
averaging of the u and v components. Hansen states that 


the computations are always stable for A, values greater 


h 


than zero. He also usesS an averaging coefficient, a, 


march as a function of AL and of the local and temporal 


step length. (When A, =0, A= Is) 
The tides, being long waves, extend their influence 


over the depth of the oceans. The equations, as stated, 
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are integrated from the sea floor to the sea surface. This 
reduces the number of variables and unknown functions by 
one. The two horizontal components of the averaged 
velocity over depth remain. 

Another assumption involved in using these equations 
is that the pressure is hydrostatic and that the changes 
in pressure are due solely to changes in surface water 
Mes tions Advection terms are considered to be small 
and are ignored when using these equations. This has 
been found to be a good assumption for long wate problems. 
These assumptions help bring equations (1),(2), and (3) 
into 2 more useable form. The problem of the geometry 
of distribution remains once these assumptions have been 
made. 

The motions of the ocean are difficult to solve 
analytically, especially if the equations are used in 
their general, non-linear form. Analytical solutions 
require considerable simplification of the boundary 
conditions and the equations can only be solved fora 
basin of regular shape and simple wind distribution. 

The Gulf of Maine is irregular in shape, non-uniform in 
depth, and exhibits varied wind patterns. To deal with 
these problems an implicit method for achieving time 
dependent solutions to equations (1),(2), and (3) was 
developed using a finite difference approach. The 


Secooulting Linlite adit terence appmemilatlons ate 
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velocity and water elevation components are given by: 
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The values of u ue V se and ein equations (4) - (6) are 


given by: 
RS LSU Gyen-1)+ url > + 
LO (man) = y KY at-i d+ Mel, ane] (10) 


ut(mon) tur Late ds 
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The half time step is given by 2. The depth in terms of 


ote and H is approximated by: 
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The external body forces composed of the wind stress and. 
barometric pressure anomoly are computed using the following 


formulas: 


ae 





(15) 


and 


Cy) 
Y* = co. t Q Po (16) 


The components of.wind stress are expressed as: 


n° 


x 2 
). AW» (Wx*+Wy ) ‘2. oe 


and 


a Ce) = i Wy (wy * 4 Wy*) Ve (18) 


maine bottom friction is assumed a quadratic funct®on of 
the velocity. The respective horizontal components are 


given by: 


(b) r 2% .r 2 ey (19) 
Tonys Retin d® 5s Lvaaivayle ( 


The partial derivatives in equations Gin (2), and 
fo) are replaced by finite difference approximations. The 
boundary values are then prescribed as known functions of 
time. The values of sea level and velocity are determined 
step by step by the difference equations starting from 
arbitrary initial values of the same quantities. 
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ITI. THE NUMERICAL MODEL 


The model developed by Walter Hansen is based on four 
major conditions to obtain solutions; (1) that® treme quae 
Smomotion must be in avsuitable form, (2) that themedrw an 
momces acting on the amea sof Shudy are given, (3) then 
boundary and initial value conditions are specified, and 
imethat the areasof investigation is known in regard to 
@eastal outline and depth distribution. ‘nese conditions 
are met in developing the Gulf of Maine model. 

The equations used are in the form of the finite 
wercrerence formulas given in the previous section. “The 
driving forces acting in the area, including the tides 
ae the open boundaries, the coriolis parameter, the 
direction and velocity of the wind field, the acceleration 
St fravity, and the coefficient of friction are Specrired 
to insure a uniquely determined solution. The open 
Doundaries are specified and the tidal constituents along 
those boundaries are known. The depth distribution and 
eoastal outline are obtained from U. S. Navy Chart Be 70c. 

The griva networkeuised anmeune Hansen model” is shown 
maeeonc GOp Ol Figure 2.) It consists of three separate 
sets of grid points: Ci) the water elevation, or Z 
mornts, at ule prid intersections, (2 )ethesuspoints 
located to the right of and midway between successive 
mepointss and (3)"the v grid points located vertically 
beneath and-mid-way between successive z points. Each 


of these points is located by the same (N,M) co-ordinate 
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designation. The N-axis is perpendicular to the first 
open boundary and the M-axis is parallel to that same 
input boundary. 


This grid set-up requires the preparation of three 


separate sets of arrays. The first set consists of symbolic 
faba for the Z grid points. The Z points are coded Oto 
land, -1 for a land-sea boundary, and 1 for a point 

located Aven the water. The separate u and v data decks 


are also coded O for land and -1 for a land-sea boundary 
Peene. The oceante points in the u and v data cards are 
coded with the charted depth, in centimeters, at the 

Pee pective points. 

The question arose early in this project as to the 
necessity of using three separate sets of data cards. 
The land and land-sea boundary areas are largely duplicated 
an- the three data decks. Since both the u and v grid 
points are located half way between succesSive Z grid 
points, another question arose as to the possibility of 
averaging sticcesSive Z grid points horizontally and 
vertically and using those values at the u and v points 
weeonout Sacrificing accuracy. 

Consequently, one set of data cards is used. The Zz 
Braid points are used, coding O for land points, -1 for 
land-sea boundaries, and charted depth in centimeters 
for the oceanic points. This alternate grid network is 
shown at the bottom of Figure 2. 


The use of only the z grid points for data make 
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several program changes necessary. The modifications 
consist of computing loops that average the values of 
successive Zz points horizontally and vertically to 
produce seuhuee of topography for the u and ~~ poumite 
These changes enter into the program in the initial data 
handling subroutine. 

Problems arose using the averaging process in this 
study to obtain data for the u and v grid points. The 
effect of averaging a boundary point, coded -1; with 
erther a land point, coded O; or an oceanic point, 
coded with charted depth in centimeters; is to produce 
a boundary point. This creates boundaries that are too 
wide. The result is that islands become larger than they 
really are and that narrow, restricted channels are 
closed. The latter effect is apparent in the Grand 
Manan Channel at the northeast corner of the area. 

The first few runs of the program ended with errors 
as the values for the tidal height in the Grand Manan 
Channel rose as high as one thousand centimeters. This 
is above built in program limits and is an unreasonable 
Meaght for the tide in this area. 

The process of correcting the boundary points 
involves additional arithmetic function statements to 
spot correct the points found to be in error. It is not 
useful to do any revamping of the Zz grid points. Live 
u and v points are corrected with direct statements 


changing those points that require alteration in value. 
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This program requires the changing of one hundred three 
values at u points and sixty-four of the values at v 
grid points. These changes enable the program to run 
Zo completion. 

Figure 2 also shows typical boundaries using Hansen's 
procedure and the revised method. In Hansen's method: 
at the top of the figure, a step type boundary is shown. 
fits boundary 12S drawn to the u and v points as ag 
blocks off the coastline. In contrast, the revised 
menwod ShowS a Step type Betndany that is drawn to the 
Z points to accomplish the same objective. The revised 
method offers less ee ae in adjusting the boundary 
to the coastline since it can be drawn to only one set 
ot -Srid points. 

Hansen's model has been tested in areas with one 
open boundary, such as Chesapeake Bay, with opposing 
epen boundaries, for the Straits of Gibraltar, and in 
the "along the coast" problem, for the area — Santa 
manrbara to San Diego, Californigsn.® This siete ee the 
Gulf of Maine uses two adjacent open boundaries and is 
an extension of the testing of the model to this new 
mode. 

The FORTRAN II program is shown in the flow chart 
ma rigure 3. —~Phe main program consists of ome program 
Sentcrol, five principal subroutines, and three plotting 
subroutines. The main program control calls the various 
subroutines when they are needed and stops the program 


when it is completed. 
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The program starts by reading 2n) all (ia ede 
in subroutine JO2, with the exeeptTon of stidalwdataean. 
the input boundary. This subroutine also computes the 
value of all compound parameters used in the program 
and finally, prints out the input data that was read in 
as a check on the accuracy of these data. This information 
is Stored in plotting Subroutine DO. 

Then subroutine JO4 prints the computed horizontal 
fields of water elevation, current magnitude and current 
feeection at the desired time intervals. These values 
are then stored in plotting subroutine DO2. 

Subroutine JO5 performs the main computations 
required in the model. It averages and smoothes the u, 
meeeand 2 values. it computes the actual u,v, and z@ 
values at each step. It calls the wind current sub- 
moucaine. The wind current subroutine, SO1, reads in 
the values of the wind fields. In the integrated 
equation wind currents are assumed to reach the bottom 
in shallow beach areas, and to extend downward to the 
top of the thermocline in regions of deep water. 

The dotted lines represent a second sequence of 
printing ’out the computed values by JO4 and the storage 
of these values for future plotting by subroutine DO2. 
The program then proceeds to compute the values for the 
special points in subroutine SO8. The program ends by 


SOS tiio tnhonia lowe ploUuline  sumpmouLtine DOG. 
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The information from the three plotting subroutines is 
then stored on a magnetic tape which in turn serves as 
the input to a Separate plotting program. 


Appendix A contains an identification of the symbols 


used in this program. 
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PV oe tHE INPOUr RDS 


The first parameter chosen is the grid si Zeneca 
Size selected depends on the area of study, the detail 
mesired, andeaiite eeu Apiee Sem mtkber core storage, memeny. 
The large area of the Gulf of Maine makes the use of a 
large array computer, the CDC 6500, necessary. The grid 
size of six nautical miles used in this study requires 
> ,000 bytes of Core memory to run the entire  prezram. 
This is over one-half the available memory of the computer. 
A smaller grid size, of say three nautical miles, requires 
even more memory. Because of rough topography a larger 
grid size results ina less accurate reproduction of the 
tides and tidal currents due to over smoothing of the 
depth distribution. 

The selection of A grid Size fixes thiewsi 7c fosmanac 
Somputation array. A Six nautical mile grid size used 
in an area five degrees:of latitude by five degrees of 
mongitude means the computational grid array will be 
51 x 51 units. The number of points on each open boundary 
is fixed by the geography along the boundary. In this 
region both open boundaries contain only points located 
over water. Each open boundary then contains 51 points. 

The selection of the special points; at which the 
effects of wind on tidal height, speed, and direction 
are studied, is based on many factors. The principal 


reason is the interest value of the point. Any of the 
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approximately 2,000 grid points Located over water are 
suitable for selection. The points chosen for this 
model are; the open sea outside Base our Massachusetts, 
the center of Cape Cod Bay, the center of Buzzards Bay 
on the western side of the Cape Cod Canal, the river 
mouth near Portsmouth, New Hampshire, the indented coast 
near Portland, Maine, the port of Rockland, Maine, Bar 
Harbor, Maine, the center of the Grand Manan Channel, a 
point in the center of Wilkinson Basin, the channel 
between Nantucket and Chappaquidick Islands, a point 
seaward of Cape Cod , a point over the continental slope, 
two points in the region:of Nantucket Shoals, and finally 
a point south of Grand Manan Island over the continental 
shelf. These points are shown in Figures 8, 9, and 10. 

Each type of topography, water depth, and coastal 
location is represented by one of these special points. 
The objective is to describe the region by comparison 
and contrast at the various Special locations. 

The acceleration of gravity is 980.665 cm/sec” fone 
this problem. This is a widely accepted value for 
gravity. 

The program provides two parameters for computing 
the wind effects on the tides and tidal currents froma 
given barometric pressure distribution. These parameters 
are the coefficient of geostrophic wind, 0.65, and the 
average air density given by 1.1627 x 107? er emcee These 


parameters do not enter into the program when uniform 
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wind fields are used. This model uses uniform wind 
fields, 

The program provides the option of using four separate 
tidal constituents at each open boundary or of using only 
the velocity of the M2 tide at these same boundaries. 

This model uses four separate tidal constituents at each 
open boundary. 


> 


The Goriolis parametem used is Siven as Gesoc7emeanies 
eo This 1S an average for the latitude of the region. 
ihe coefficient of sboetrom friction, the propemtcma lama, 
factor between the wind stress and square of the wind 
velocity, and the smoothing parameter used with the 
Moet rticient of horizontal eddy diffusivity are piven by 
values which Hansen found to be the most suitable after 
much experimentation with the model for widely separate 
Meates Of water. The coefficient of friction used an 
this study is 0.003. This as a dimensionless parameter. 
The components of the wind stress and the square of the 
wind velocity are related by the proportionality factor 
Die xX 107, The numerical smoothing parameter, a, is 
0.998. Hansen obtained the most accurate results with 
his model using these values. 
The one-half space step is determined by the grid 
rzec.ee in thissproblem it itowedven a@se5. 5500 x 10° 


centimeters, or three nautical miles. The maximum length 


of the time step is determined from the maximum grid size 
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and depth in the area according to the Courant-Friedrich- 


Levy criterion. 


42 24 


sl Cq Hmax> 


ine thas formula, t is the (tame ain seconds , the grid 





My | (20) 


Size in centimeters, g the acceleration of gravity, and 
a the maximum depth in the area of computation in 
centimeters. When this criterion is applied to the Guif 
of Maine, a time step of 31 seconds results. 

The wind field characteristics are determined by 
the -programmer. The wind fields are specified by the 
wind velocity in m/sec, the wind direction computation 
co-ordinates, and the time, in seconds, when the wind 
starts Bowne: The wind fields constitute one of the 
two primary inputs to the program. Their values are 
introduced at each grid point for every time step 
throughout the program. 

The tidal constituents are obtained from tabulated 
values for various locations throughout the world. The 
M2, K1, S2, and N2 constituents for Grand Manan Island 
are used for the eastern open boundary and the M2, K1, 
52, and N2 constituents for the western end of the Cape 
Cod Canal are used for the southern open boundary. 
finese values are introduced at each respective open 
Poundary point for every time step throughout the program. 
Table I summarizes the values of the tidal constituents 


used. 


DD) 





SrTnrTTETTE ESET nnn 
TABEY a 


a 


Grand Manan Island -W. Cape Cod Canal 


Gonstituent Speed Direction Speed Direction 
(cm/sec ) oa (cm/sec ) a 

Be EE ——————eeeee— 
M2 2 aS 339 510) 6) 230 
S2 BOr@ 015 10.0 20S 
K1 omy 131 10510 O80 
N2 EOL Bae 15.0 205 


ee eres 


The final data decisions are those concerning how 


many runs to make and the various time considerations 


involved with each run. A total of seven separate runs 
are reported in this thesis. The first is a run using 
only the tidal inputs and no wind. This run requires 


a total of approximately five hours of computer time on 
the CDC 6500 to accomplish thirty-five hours of computations. 
Hnais is sufficient to cover almost three semi-diuznal 
tidal cycles and to establish the tidal circulation. 

The next four data runs use 20 knot wind fields 
m@rom the four cardinal directions. In each case the wind 
is introduced at the five hour mark and terminated at 
thirty-five hours. Output from the model is obtained 
at two hour intervals. The objective is to combine the 
wind and tidal circulation and see how the wind fields 
modify the tides and tidal currents. 

The final two data runs use the same format for 


length of computation, introduction of the wind, and 





termination of the wind. The velocity is increased To 

30 knots and the directions are altered from the southeast 
and northeast respectively. The input parameters and 
their data card format and arrangement are contained in 


Appendix B 
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V. RESULTS 


The results of this study consist of a computer sour 
put for each of the seven data runs. Each computer print- 
out contains the water ele aetaion in centimeters, the 
resultant tidal current speed in centimeters per second, 
and the tidal current direction in degrees true for each 
Srid point at two. hour intervals. Each run contains 
nineteen such printouts. The same information is printed 
aut for each of the fifteen special points. 

This produces a great deal of data. Each two hourly 
printout of values contains computed information for nearly 
two thousand oceanic grid points. It is difficult to 
interpret all of this information without some reasonably 
detailed plots. The ideal situation is to develope a 
computer program to plot part or all of the data manually 
in some suitable format. Because of difficulty with the 
peot program and time limitations the alternative of 
plotting part of the data is necessary. 

The results are discussed in two sections. The first 
section concerns the verification of the Hansen model in 
the Gulf of Maine. The second section discusses the 
effects of various wind speeds and directions on the 
tidal water level, tidal current speed, and the tidal 


current direction. 
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A. VERIFICATION OF THE HANSEN MODEL 

The verification of the model requires that some 
observed data be compared with the numerical values 
obtained from the model. The observed data used in this 
mostance is tidal current data compiled by the U.S. 

Weast and Geodedic survey. is data represents current 
measurements extracted from light ship observations over 
eapperiod of several years. The currents are not Correlated 
with the winds present at the times that they were 
measured. This means that the comparison of data that 
follows actually compares average currents produced over 

a period of years by widely varying atmospheric conditions 
with the instantaneous values produced by the model for 

a Unique set of prevailing conditions. 

Figure 4 presents tidal current speed in knots and 
ferdal current direction aS indicated by the vectors, for 
the observed Coast and Geodedic Survey data. This 
description pertains to the time of five minutes after 
low water at Boston, Massachusetts. Immediately beneath 
this data, Figure 5 presents the same information 
obtained from the model for the same approximate locations 
and ae of low water at Boston, Massachusetts. In 
extracting comparison information from the model, the 
closest grid point to the location of the observed data 
is used. The measurements are not compared at exactly the 
Same geographic location. Figures 6 and 7 perform the 


Same respective functions for the high tide situation. 
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Comparison of these figures shows broad subjective 
aereenenl ssetuweene lire Ouse wed data and that computed by 
Hansen's model for the high and low tide situations. The 
agreement oe is general and not in detail. “in fact, 
there are specific points at high and at low tide where 
disagreement is evident. The tidal speeds do show an 
order of magnitude agreement. The directions of the tidal 
currents appear reasonable for both high and low tide. 
The agreement shown suggests that currents in the region 
are probably tidally dominated. 

The information used in Figures 5 and 7 is obtained 
from the no wind run of the model. The reference station 
used is Boston, Massachusetts to conform with the Coast 
and Geodedic Survey data. 

The region of comparison is known as George's Bank 
and covers an area bounded by latitude 40°30'N and 42°N 
and mene ee 67°30'W and 7O°W longitude. This represents 
only a portion of the entire Gulf of Maine area. 

The meteorological conditions prevailing in the Gulf 
of Maine vary considerably from month to month. There 
Mano ticeable daidperence: ime ther pattern of the wmean 
wind flow between winter and summer. In February, a 
month typical of the winter regime, the mean winds are 
westerly with a velocity ranging from 7 to 27 knots. The 
Summer season, represented by August, has mean winds that 
are Cio metre by southerly to southeasterly with 


velocities from 3 to 16 knots. The Gulf of Maine isa 
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reelon, of cyclonge storm a cede and exhibits a wide 
range of wind speeds and directions during any given 
mer iod, 

The verification if Somewhat disapporm ten eee ae 
many questions remain unresolved. The principal obstacle 
remains ascertaining the conditions prevailing when the 
observed data was recorded by the Coast and Geodedic 
Survey. Unfortunately, the data does not provide the 
meteorological conditions prevailing at the time the 
measurements were made. Since nee are averages obtained 
over many years it is assumed that many separate wind 
fields prevailed during these current measurements. This 
appears to preclude a more accurate verification without 
additional specific information. 

Figure 8 presents the tidal current circulation in 
the Gulf of Maine when high tide occurs at Boston, 
Massachusetts. Figures 9 and 10 show the same information 
for slack water and low tide respectively. Slack water 
and a zero tide level occur almost simultaneously for 
this location. All three diagrams are compiled from the 
no wind data run and show tidal circulation not driven 
by any specific wind. 
roe LH EFFECT OF WIND FIELDS ON TIDAL HEIGHT, TIDAL 
CURRENT SPEED AND TIDAL CURRENT DIRECTION 

die ettec, of the var 10us winedetselds on the tides 


mia tidal weurirent Cauculation as prasented in a Series of 
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tables. The values used are those Computed for thewtatteen 
special points for each of the seven data runs. 

The duration of the wind is considered by choosing 
times such that the fully arisen sea and non-fully arisen 
sea cases are considered. The wind is introduced at the 
five hour point in each data run. The first time chosen 
is at 10 hours or five hours after the wind starts. Five 
hours is not sufficient time for a twenty or thirty knot 
Totcdeto build inte a fully arisen sea (i.e. time for the 
energy content of the sea to be saturated). The second 
elected time is at 24 hots Gr 19 hours after the ounces 
of the wind. This time frame is sufficient for the 20 
knot winds to generate a fully arisen sea but not for the 
30 knot winds to accomplish the same result. The final 
time selected is at 35 hours or 30 hours after the onset 
of the wind. This allows the 30 knot winds to generate 
aQeprully arisen sea and explores the effects of allowing 
the wind to blow over an already fully arisen sea for 
Pmeoven hours in the case of the 20 knot wind fields. 

ane wager welevation atten 10 hours, (5 hotsss som wind) 
1S summarized for the fifteen special points in Table Il. 
Simalar information for 24 hours and for 35 hours of wind 
is contained in Tables III and IV respectively. Table V 
contains the maximum change in water height produced by 
any wind field at each point for times ten, twenty-four, 
.and thirty-five hours. The range of maximum wind produced 


changes runs from no change at New Bedford, Massachusetts 
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TABLE V 









MAXIMUM CHANGE 
IN WATER ELEVATION (CM) 


24 HOURS 2D HOURS 
BOSTON, MASS. es aT 


NEW BEDFORD ,MASS. ; lesen 


PORTSMOUTH ,N.H. (omcs 11.48 
PORTLAND, ME. iso et 9.35 


ROCKLAND ,ME. 2.92 


BAR HARBOR,ME. 
GRAND MANAN IS. 





LOCATION 
ie Se Wres 





HERE : 


CASHES LEDGE Be ae 
CHAPPAQUIDICK IS. ey eee 
OE ee | 
CONTINENTAL SLOPE ie leno | ioe 


NANTUCKET SHOALS Di Oe 


REAT SOUTH eSsHOALS 


R. MANAN CHANNEL 
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and Grand Manan Island to a change of 18.4Gecentamete:s> 

for Chappaquidiek Island at the 2% hour time. Mii hemeneatest 
change in tidal height due to any ean field amounts to 
78.46 centimeters Or approximately one haltila Gootemsounenme 
locations and conditions specified in this study. In 
computing the differences reported in Table V, the no 

wind data ivom Tables Ti ttiweand IV serves as tie seaec 

or reference oe No attempt is made to identify the 
wind field responsible for producing the greatest change 

at each Pocite 

pWe= changes at deepewater slocations are ysmal ?mxliske 
Ane eréatest changes in tidal height induced by the wind 
occur at points in relatively shallow water. This may be 
a consequence of the averaging of the model equations 
over depth. The lack of any change at New Bedford, 
Massachusetts is attributed to the effects of geography - 
mars location iS Sheltered by Doundary points from every 
wind field used in this study. The maximum changes 
observed do pot necessarily occur when the wind has been 
peeowlne for its full duration. 

The effects of the various wind fields on the tidal 
current speed and direction at the special points are 
presented in Tables VI, VII, and VIII for the ten, twenty- 
four, and thirty-five hour times. <A summary of the 
maximum changes in these quantities is presented in Table 


- TX. Once again only the greatest changes are considered 
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TABLE Ix 


MAXIMUM CHANGE 

IN 
TIDAL CURRENT MAGNITUDE (CM/SEC) 
TIDAL CURRENT DIRECTION (°T) 


10 HOURS | 24 HOURS 35 HOURS 
BOSTON, MASS. 49/53 46/165 
CAPE COD BAY 41/174 
EW BEDFORD ,MASS. O07 176 
RAND MANAN IS. 0/0 0/0 0/0 — 
CONTINENTAL SLOPE 13/45 ye 


2eInG 11/125 












BAR HARBOR,ME. 





NANTUCKET SHOALS 







GREAT SOUTH SHOALS 7/9 9/13 Lig 36 
R. MANAN CHANNEL 2s of) 
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and no attempt is made to identify the particular combination 
of wind speed and direction that causes the greatest change 
at a given location for a given time. 

The variation in tidal current speed runs from a high 
mt oO cm/sec at Boston, Massachusetts at 24 Nourcmccme 
low value of only 1 cm/sec in the Grand Manan Channel at 
the 10 hour point. Thus the range of changes varies from 
imooeot a knot to 1.2 Knots in speed. 

The variation in the direction of the tidal currents 

due to the addition of the winds is great. The range is 
froma oe small change of only 2 degrees to an almost a 


complete reversal of direction of 178 degrees. 


TABLE X 


WIND VELOCITY (KNOTS ) 10 20 30 4O 50 


AVERAGE CURRENT VELOCITY DUE 

fOmWwiIND AT FOLLOWING LIGHT- 

SHIP STATIONS 
BOSTON AND BARNEGAT OFT (gO SeeOP2 a0 eeons 
DIAMOND SITOAL (CAPE HATTPRAS )»O.5 0.6. “Omgeen Orc mao 


ALL OTHER LOCATIONS OZ O<3 Out 0.5 0.6 


fatal eeurrent tables published by the Us )> .Coacy and 
Geodedic Survey discuss the effect of the wind on tidal 
Strona velocity.. Pable xX ShowsSrtime 1epor Ged souseryved 
differences caused by various wind speeds. These data 
are based on observations made by lightships during a 


period of 54 years. They do not consider the duration of 
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the wind field and this helps explain the differences 
between the data reported in Table IX and the observed 
values shown in this table. 

Points located over deep water, such as the Cashes 
Ledge point, exhibit very little change in tidal current 
speed or direction due to the effects of wind. The 
largest response of these quantities to different wind 
fields eneeaee atepoints located in Shallow watverwarcace 
Geographic position again enters into the evaluation of 
these factors. The largest changes in speed ain direction 
of the tidal currents occur at points near the boundary 
between the land and sea. Such points as Cape Cod Bay, 
New Bedford, Massachusetts, Rockland, Maine, and Chappaquidick 
Island all exhibit large deviations in the direction of the 
tidal current in response to various wind fields. 

The use of equations integrated over depth spreads 
out the effect of changes at deep water locations over 
the full depth at those points. All equal amount of 
change at a point in more shallow water is distributed 
over a lesser depth. This helps explain the greater 
changes noted in shallow water locations. 

The model also considers that there is no net mass 
transport perpendicular to the shoreline. This introduces 
the explanation that tidal currents flowing at the surface 
are offset by a subsurface or bottom flow in the opposite 
direction, so that the net mean current as presented by 


this model 1s Zero. 
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VE. CONCLUSIONS 


The conclusions reported in shale study are divided 
mato three major categories. The firs comecrns siie ase 
of Hansen's model in a mode with two adjacent open 
boundaries and a modified method for handling the topographic 
input data. The next category is the verification of the 
marasproduced by Hansen's model with observed data.) Mec 
final category concerns the effects of wind fields on the 
tidal heights, tidal speeds, and tidal current directions 
in the Gulf of Maine. 

The use of Hansen's model with two adjacent open 
boundaries produces reasonable results for the Gulf of 
Potie region. The modified grid network used in this study 
requires minor adjustments in its use but is workable and 
fmesoults in less time and effort expesided in coding and 
key punching topographic data. 

Hansen's model may be used in a variety of ways. The 
fact that the use of the model with adjacent open boundaries 
produces reasonable results does not mean that these are 
the best answers possible. The model can be used with one 
open boundary, various parameters are subject to change, 
and any combination of wind speed and direction is possible. 
The verification of the Hansen model reported in this study 
is only in broad subjective agreement with observed data. 
Poe verification 1S for a Pivenmset of Conditions whieh 


can be altered to produce an entirely different result. 
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Lt is Surfiertent sve say jthac eee the Hansen model with 
adjacent open boundaries produces reasonable results 
which are in general agreement with some observed data 
for the Gulf of Maine. 

The conclusions concerning the effects of various wind 
fields on tidal height, tidal current speed, and tidal 
Current direction are limited an ©s<cope and are wmeccsscarai., 
Goniined to the fifteen special points used in the study. 
Any attempt to extend these conclusions to the other 
points in the Gulf of Maine is eee aren as each point 
must be looked at individually. 

The effect of the wind fields used in this study on 
the tidal height at the points in question is slight with 
me lareseste change being only one half a foot. Ihe 
greatest changes are found at points located in shallow 
water. In contrast the smallest changes are found for 
points located in deep water. Geographic location appears 
non ibe important in considering the effect of wind fields 
oo the height of the tides. | 

The efiect of the wind fields on tidal current 
velocity at the specified locations used in this thesis 
varies considerably. The wind speed and direction is 
found to have less effect at deep water locations than it 
does at shallow water points. This is expected because 
the equations used in this study are integrated over depth 
and hence, the effects of the wind are spread out over a 


large area in deep water, vice a much smaller area in 
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shallow water. Physically this is a reasonable result. 

The speed of the tidal current increased by 1100% amounting 
to 1.2 knots at Boston, Massachusetts for the greatest 
ehange recorded. This may be of sigmiiicance Vana sindicates 
that the wind speed and Pee does influence tidal 
current speed. 

The effect of wind speed and direction on the resulting 
pe ccit ion of the tidal currents is pronounced. The 
variation is again found to be greater for points located 
in shallow water than for deep water locations. In some 
cases, the direction of the tidal current is found to 
almost completely reverse under the influence of a wind 
field. This indicates that wind speed and direction can 
not be ignored in the prediction of tidal current direction. 

LG ee aeeneesne to realize that the current experienced 
by an observer at a given point is a summation of the 
mudal current, wind driven current and any major oceanic 
eurrent, such as the Gulf Stream, at that point. This 
study concerns itself entirely with tidal currents and 
mire eCtfects of various wind fields acting on them-: 

The Hansen model provides a reasonably fast, generally 
mechiabve MmeLMod of predicting tidal current speeduand 
Pereco1on, and tidal heient for any fiven Unrquely specified 
Sond’ tLOns ae In Unas refrard, “there 1s a -preag neeq for 
faba wi tie Wwi Hen bemwer Diy Sthne model wn CUur remesiieasunt Cmcnita 
should be combined with the atmospheric conditions pre- 


vailing at the time they are made to enable investigators 
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to study the complete atmospheric oceanic system of 


interactions. It appears the model is most useful for 


storm surge prediction. This is vital information for 


coastal engineers and residents of coastal areas. The 


model is capable of producing information for any point 


specified 
knowledge 
Emeced, Or 

Laas 


mie tides 


and hence can assist an falttinge the vord vor 

for areas where no observed data on tidal height, 
directions 1s taecorded:. 

project represents only one attempt to study 


and tidal currents in the Gulf of Maine. As 


such it only begins to explore a most complex and involved 


subject. 
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Vii. SUGGESTIONS FORSEULURE AESeaker 


The Hansen model is capable of many separate modes 


of operation. Use of the model in the Gulf of Maine region 
should be continued and the model should be used in as 

many different modes of operation as feasible in order 

to find by comparison and contrast that combination of 
boundary conditions and parameter oes that yields 

the most realistic results for the region. 

The rhea: grid data handiing method requires 
Werifacation. The Hansen model should be run again for 
the Gulf of Maine with the same parameters using the 
original method of handling the topographic data. A 
comparison of results between the two methods should 
help determine the worth of the modified method of 
handling these data. 

The rugged, fractured topography of the area suggests 
fMmiak a £rid Size of Six Nautical miles may be too course. 
A grid size of three or maybe even two nautical miles may 
be necessary to bring out more accurate results for the 
tidal height, tidal current speeds, and tidal current 
directions in this area. The smaller grid size should 
elamainate any efte€ets of over smoothing of the depth 
eavstribution possibly cavlced by the six nautical mile 
grid size. This may require the sub-division of the Gulf 
of Maine into two or more Seen regions in order to 


reduce the computer time involved in such efforts. 
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The investigation of the effects of variocus wind 
fields on the tidal heights and tidal currents is worthy 
of further efforts. This is essentially storm surge 
prediction and it is of practical interest to mariners, 
coastal engineers, and many other groups of persons. It 
could, some day, prevent property damage or even loss of 
life by adequately predicting the effects of a severe 
storm. 

Constant experimentation with the Hansen model is 
the only way to discover whether it is in fact valid for 
all types of water regions. it's use should be extended 
to Es Many SeCparatewoceanic areas as Leasiole. 

The Hansen model used in this study is limited in 
the sense that it does not accurately portray the density 
distribution of the ocean. The recent development and 
testing of a so-called multi-layer model should help 
overcome this difficulty. The results of the multi- 
layer model and this single density Hansen model can 
someday be compared to determine which most accurately 
describes the tides and tidal currents in the Gulf of 
Maine. It may be that for some regions the single density 
model is sufficiently accurate. in contrast, the mulva— 
layer models may offer significant improvements in the 
eccuracy Of Teproduicing the tides and tidal currents. 
Much remains to be accomplished before mother nature is 


duplicated by a deck of computer cards. 
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APPENDIX A 


ABBREVIATIONS AND PARAMETERS USED IN THE PROGRAM 


M 

N 

Z(N ,M) 
U(N,M) 
V(N,M) 
RAD(N,M) 


ANG(N,M) 


QU(N,M) 


QV(N,M) 
REST(N,M) 


DIR(N,M) 


ZM(N,M) . 
ZST(N,M) 
HTZ(N,M) 
HTU(N,M) 
HTV(N,M) 
HGU(N,M) 
HGV(N,M) 


XK(N,M) 
YK(N,M) 


A(T) 


id index para lle Wewomemigeecumices) 


grid index (perpendicular to entrance) 
water elevation (CM) 

U-component of velocity (CM/SEC) 
V-component of velocity (CM/SEC ) 
resultant current speed (CM/SEC) 


angle from geographic north of resultant 
current speed 


fields for summation of the U and V components 
of currents for computation of the rest 
currents 

see QU(N,M) 

average speed of rest current (CM/SEC) 


direcuilon Otis es tweurtr ent (geographic co- 
ordinates ) 


used in JO5 as a smoothing intermediate 

used in JO5 as an averaging intermediate 

water depth at the water elevation (Z) points 
water depth at the U-points 

water depth at the V-points 

actual water depth at the U-points (HGU=HTV+Z) 
actual water depth at the V-points (HGV=HTV+Z) 


U-component of the wind current divided by 
depth 


V-component of the wind current divided by 
depth 


Chawactue: 1stics of the wind field 
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A(1) 
alee) 
A(3) 


NU(T) 


MU(T) 


NZ (1 ) 


MZ(1) 


NK(1) 


MK(T) 


NA(T) 


MA(T) 


ZA ( 1) 


mae 1) 
U1(T) 
we (T ) 


73( I ) 


ZA T ) 
U3(T) 
mer ( 7 ) 
Vi(1) 
V2(1) 


V3(TI) 


time when the wind starts (SEC) 

wind speed (M/SEC) 

wind direction (computation co-ordinates) 

wind field delimeters (N and M co-ordinates 
of the upper and lower right corners of 
the wind field) 

see NU(T) 


co-ordinates (N,M) of the points used as 
selected special points 


see NZ(I1) 


co-ordinates (N,M) of the points on the first 
open boundary 


see NK(1) 


co-ordinates (N,M) of the points on the second 
open boundary 


see NA(T) 


amplitudes of four tidal constituents at the 
first open boundary (CM) 


see Z1(T) 
see Z1(1) 
see Z1(T1) 


amplitudes of four tidal constituents at the 
second open boundary (CM) 


see Z3(1) 

see Z23(1) 

see 73(1) 

names of special points 

water elevation at special points (CM) 


speed of the current at the special points 
(CM/SEC ) 
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V4(T) 


AFGN 


> TGMA 
ALPHA 
R 

ROL 
RBETA 
C 

DL 
DT 

al 

TI 

2 

NE 
ME 
TE 
KKE 


KO 


IZE 


NG 
LUE 
TW 
(hake 


LEN 


direction of the current at the special points 
arbitrary problem number 

coriolis parameter (1 /SEC) 

acceleration of gravity (CM/SEC* ) 

angular velocity of the M2 tide (RADIANS/SEC) 
smoothing parameter 

friction coefficient 

air density (GM/cM? ) 

coefficient of geostrophic wind 

drag coefficient 

+ step in space (CM) 

> step in time (SEC) 

time (SEC) 

interval between printouts (SEC) 

field cut puumeaunter 4 SEC) 

field size delimeters 

see NE 

end time of computation (SEC) 

number of wind fields 


the number of points on the open boundary 
minus one 


number of open boundary points (first open 
boundary ) 


number of points on second open boundary 
twice the number of wind fields 

time when the wind starts 

counter 


see TIC 
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POC 
BETA 
Al 
A2 
A3 
A4 
A5 
C1 


C3 


SL 
NURU 
NURV 


JA 


NEH 
MEH 
N EHH 
M EHH 


WERTO 


WERTU 
WERTL 


WERTR 


see TIC 
(1-ALPHA ) /4 


2DT 


F(A1) 


R(A1 ) 

De, Die 

G(A4) 

Of Sivoo wala) ore ew tcl 

Cc(A(2))x(10,000 is a unit conversion factor 
which converts (M/SEC) to (CM/SEC) after 
A(2) is squared 

time when wind stops 

number of ‘special points 


printout line counter 


indicator, sac .OMiset Us —1\=— 7 =O rme ine Hn 
read initial values 


NE-1 
ME-1 
NE-2 
ME-2 


storage location for U and V during smoothing 
and averaging 


see WERTO 
see WERTO 


see WERTO 


WERTOL 


WERTOR 


WERTUL 


WERTUR 


sec 


=o Cc 


Ssce 


see 


WERTO 


WERTO 


WERTO 


WERTO 
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WURZEL 


GRZ 


SQRT(ZM*-ZST*) 


A3(WURZEL ) 
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APPENDIX B 


INPUT PARAMETERS AND THEIR FORMATS 


pie wData InpuleCards 


Garda 1 


Fommat 24ie 


JA,NE,ME,1IZE, IUE, KKE, NURU,NG 


JA 


NE,ME 
IAD 
IUE 
KKE 
NURU 


NG 


ward 2 


indicator. JA=Crme tees) =) Vee Os A= leer ead 
initial values 


field size delimeters (51,51) 

number of points at first open boundary (50) 
double the number of wind fields (2) 

number of wind field characteristics (3) 
number of special points (15) 


number of points at the second open boundary (50) 


Format 9RGa3 


AFGN,G,ALPHA,RBETA,C1 


AFGN 
G 
ALPHA 
RBETA 


C1 


Card 3 


problem number (7001) 

acceleration of gravity (978CM/SEC“) 
smoothing parameter (0.998) 
coefficient of geostropic wind (0.65) 


Wild In@tearom-wO-moswind., (=wing 


Format 9F8.O 


Mie Tee Ty Dinc2, Sa 2,73 


DT 


TE 


TW 


T1 


= time step (15 SECS) 
length of computation (126000 SECS) 
time when wind starts (18000 SECS) 


interval between printouts (7200 SECS) 


7O 





Card (cont. ) 


a2 


SI 
T 


T3 


Card 4 


field output iccumter (0 if outputs are desired 
from the start of the problem; otherwise 
any other delayed starting time, e.g. 7200 
SECS ) 

time when the wind stops (126000 SECS) 


time (initialized at zero) 


time when plots are desired (every L200 SECS ) 


Format sob ile 


DL,F,SIGMA,R,ROL,C 


DL 


F 


SIGMA 


R 


ROL 


ig@ard 5 


NZ MZ 
NZ 


MZ 


Card 6 
NK ,MK 
NK 


MK 


ea rar 7 


NU , MU 


NU 


half the grid size (CM) 
a 


Cordolasepanametera(oes > clon 
angular velocity of the M2 tide (1.4088 x 107") 
fPpict ton eoe hhc emt | 02003) 

denis t ys oueitinemecna len dG2 7x 1077) 

°) 


drag coefficient (3-2 x 10° 
Pormat 23 

N coordinate of special points 

M coordinate of special points 

Format 2413 

N coordinate of points on first open boundary 
M coordinate of points on first open boundary 


Format 2413 


N coordinate of wind field delimeter 
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engel (ecrocaitim) 


MU 


Card 8 
NA ,MA 


NA 


MA 


Card 9 
v(T) 


Bel) 


Card 10 
rr } 
A(1) 
A(2) 


AND) 


ard 11 
HTZ 


HTZ 


Card 12 
Cae Ng Ue 
Z 1 

Z2 

U1 


we 


M coordinate of wind field delimeter 


Format 2413 


N coordinates of points on second open 
boundary 


M coordinates of points on second open 
boundary 


Format 1045 


names of selected special points 


Pormat Jha. 


time when the wind starts (18000 SECS) 
wind speed (M/SEC) 
wind direction (in degrees in computation 


coordinates ) 


Pormat. 1266.0 


depth cards (usually more than 10 cards; 217 
in this program) 


Format 9F8&.3 
amplitude of first component (CM) 
amplitude of second component (CM) 
ampli tcude Tot -thard componens: (CM ) 


amplitude of fourth component (CM ) 


Ge 





Gan dans 

1 fae ry Aare is en OL 
Z 3 

Z4 

U3 


U4 


Format 9F8.3 


amplitude 
amplitude 
amplitude 


amplitude 


of 


of 


of 


of 


first component (CM) 


second component (CM) 


third component (CM) 


fourth component (CM) 
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